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ABSTRACT 

We use cosmological simulations from the Aquarius Project to study the orbital history of the 
Large Magellanic Cloud (LMC) and its potential association with other satellites of the MiUcy 
Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the 
LMC position and velocity at either of its two most recent pericentric passages. This suggests 
that the LMC is not necessarily on its first approach to the MW, provided that the virial mass 
of the Milky Way is as high as that of the parent Aquarius halo; M200 = 1.8 x 10^^ Mq. The 
simulation results yield specific predictions for the position and velocity of systems associated 
with the LMC prior to infall. If on first approach, most should lie close to the LMC because 
the Galactic tidal field has not yet had enough time to disperse them. If on second approach, 
the list of potential associates increases substantially, because of the greater sky footprint 
and velocity range of LMC-associated debris. Interestingly, our analysis rules out an LMC 
association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell 
& Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that 
the direction of the orbital angular momentum is a powerful test of LMC association. This test, 
however, requires precise proper motions, which are unavailable for most MW satellites. Of 
the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly 
associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and 
Sculptor rule them out as potential associates, but this conclusion should be revisited when 
better data become available. The dearth of satellites clearly associated with the Clouds might 
be solved by wide-field imaging surveys that target its surroundings, a region that may prove 
a fertile hunting ground for faint, previously unnoticed MW satellites. 

Key words: galaxies: haloes - galaxies: formation - galaxies: evolution - galaxies: kinematics 
and dynamics. 



1 INTRODUCTION 

The Large and Small Magellanic Clouds (LMC and SMC, re- 
spectively) are unusual satellite galaxies. They are exceptionally 
bright, and so close to the Milky Way (MW) that recent studies 
have concluded that fewer than one in ten MW-like systems are 
expe cted to host satellites with properties similar to the Clouds 
jBovlan-Kolchin et al, 2011; Btisha et al. 2010; Lares et al. 2oTl|; 
iLiu et aljlioni : Iguo et allBoi ll : iTollerud et al.ll201li) . The short 
crossing time at their present Galactocentric distance suggests that 
they may have already completed a number of orbits in the Galac- 
tic p otential (Murai & Fujimoto 1980: Lin & Lvnden-B"eiilll982l : 
iGardiner et al. 1994 : van der M arel et al. 2002i) while their prox- 
imity suggests that they form a bound pair. 

If the LMC/SMC are truly physically associated, then it is 



likely that they were once pa rt of a larger system : the "Greater 
Magellanic Galaxy", to quote Lvnden-Belll ( 119820 . This idea has 
prompted searches for evidence that other satellites might have 
been in the past associated with the Clouds. This is encouraged, in 
part, by the "polar" distribution of the brightest Gala ctic satellites, 
which seems to trace the orbital path of t he Clouds dLvnden-Belll 
ll976l) . lLvnden-Bell & Lvnden-Belll ( Il995t) . for example, suggested 
a possible association between the LMC, SMC, Draco, Ursa Minor, 
Carina and Sculptor as part of a common "Ghostly Stream". 

The association between different satellites has recently re- 
ceived renewed attention, motivated mainly by coherence in the 
pos ition and velocities of satellite pairs such as Leo IV and Leo 
V ( lBelokurovet"aLl 12008*)) and of satellites near the Sagittarius 
stream; e.g.. Segu e 1 (Niederste-Ostholt et al. 2009), Bootes II 
jKoch et all [2009I) . and Segue 2 jBelokurov et alj |2(X)9|) . These 
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ideas have received support from tlie realization that most satel- 
lites should have been accreted into the Milky Way as part of 
multiple systems and have gained momentum because testing the 
predictions of such scenari o has become pos s ible using realistic 
cosmological simulation s i Sales et al] 2007 al : Li & Helmil 20081: 
D'Onghia & Lake! l2008l : iLudlow et alli2009l : iKlimentowski et all 



20101) 



The Magellanic Clouds are also unusual in their kinematics. 
The latest pro per motion measurements of stars in the La rge Mag- 
ellanic Cloud jKallivavalil et alJl200d : rPiatek et alj|2008l) indicate 
that the LMC has a much higher tangential velocity than previ- 
ously thought (Vf ~ 370 km/s), raising questions as to whether 
the LMC and SMC ar e actually bound to each other or even to 
the Galaxy as a whole. iBesla et al J ( l2007h . for example, have used 
the new kinematic data t o revise earlier orbital models (see, e.g., 
Ivan der Marel et 31.1 12002*) and concluded that the LMC and SMC 
must be on the first pericentric passage of their orbit around the 
Milky Wa y if the MW virial m ass is of the order of ~ 10^^ Mq, as 
argued bv lKlvpin et alj ( l2002h . 

The discussion above depends sensitively on the assumed 
virial mass of the Galaxy. Bright satellites are more common 
around more massive primaries (Guo et al. 2011), and higher pri- 
mary masses make it easier to accommodate high orbital speeds for 
the satellites. Even a facto r of two increase in virial mass can make 
a difference jPiatek et al.i 2008: Shattow & Loeb 2009), which is 
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The escape speed at r = 50 kpc from a Navarro-Frenk- White 
halo with virial masfl K'hm = 2 x 10^^ Mq is of order ~ 500 
km/s, which would mean that the Cloud s are safely bound to the 
Galaxy despite their high speed. Further, iBesla et ah I (120071) show 
that virial masses as high as that imply a radial period of just about 
3 Gyr and hence the possibility that the LMC and SMC have com- 
pleted multiple orbits around the Galaxy. A first infall scenario is, 
on the other hand, compelling for a number of reasons, including 
the recent successful modeling of the Magellanic Stream as a tidal 
relic of a recent interaction betw een the clouds before infall into 
the Milky Way jBesla et alJ|20rQ) . 

The discussion above highlights the fact that basic issues such 
as whether the LMC and SMC are on their first approach, or bound 
to each other, or truly associated with other MW satellites, remain 
unresolved. These questions are clearly interrelated. For example, 
the relative positions and velocities of satellites of past LMC as- 
sociates would be very different if the Clouds are on first or sec- 
ond approach. Careful modeling is therefore required to interpret 
current data, especially because the orbit of the Clouds should be 
eroded quickly by dynamical friction and by ti dal mass los s , limit - 
ing the applicability of models such as that of iBesla et all j2007l) . 
which assume that the Clouds evolve in a rigid Galactic potential. 

We address these issues here by using the Aquarius Project, 
a series of cosmological simulations of the fo rmation of a Milky 
Way-sized halo in the LCDM paradigm ISpringel et al.ll2008l) . Be- 
cause of their extraordinary numerical resolution, we are able to 
trace the orbits not only of massive subhalos, but also of subhalos 



^ The virial mass, Af2oo , is defined as the mass contained within r2oo > the 
radius of a sphere of mean density 200 times the critical density for closure, 
Pcrit = /SttG. This choice defines implicitly the virial radius of the 
halo, r200 , and its virial velocity, V200 • 



within subhalos, thus enabling a realistic assessment of their orbital 
paths before, during, and after their accretion into the Milky Way 
halo. 

This paper is organized as follows. In Sec.|2|we provide a brief 
description of the Aquarius simulations. We analyze the orbit of an 
LMC candidate in Sec. 13. II and the spatial and kinematic properties 
of their associated subhalos in Sec. 13.21 We use these results to 
explore in Sec. 13.31 which satellites of the Milky Way might have 
been associated with the Clouds in the past. Sec.|4|summarizes our 
main conclusions. 



2 THE NUMERICAL SIMULATIONS 



2.1 The Aquarius Project 

The Aquarius Project dSpringel et aLllioOsI) is a series of cosmo- 
logical simulations of the formation of six dark matter halos with 
mass consistent with that expected for the halo of the Milky Way. 
The simulations assume the ACDM co smology, with para meters 
chosen to match the WMAP 1-year data dSpergel et aLll2003l) : mat- 
ter density parameter, JIm = 0.25; cosmological constant term, 
JIa = 0.75; power spectrum normalization, erg = 0.9; spectral 
slope. Us — 1; and Hubble parameter, h = 0.73. 

The halos were identified in a large N-body simulation of a 
cube 137 (100 h^^) Mpc comoving on a si de, a lower resolution 
version of the Millennium-II Simulation teovlan-Kolchin et al.l 
2009). This volume was resimulated using the same power spec- 
trum and phases of the original simulation, but with additional 
high-frequency waves added to regions encompassing the initial 
Lagrangian volume of each halo. The high-resolution region was 
populated with low-mass particles and the rest of the volume with 
particles of higher mass (Power et al. 2003). 

The six Aquarius halos are labelled 'Aq-A" through "Aq-F". 
Each was resimulated at different resolutions in order to assess nu- 
merical convergence. A suffix, 1 to 5, identifies the resolution level, 
with level 1 denoting the highest resolution. Between levels 1 and 5, 
the particle mass ranges from nip — 2x 10'^ Mq to = 3 x 10^ AIq. 
At z = 0, the six halos have similar "virial" mass, roughly between 
1 and 2 x 10^^ M p . For further details of th e Aquarius Project, w e 
refer the reader to lSpringel et al.l(l2008l) and lNavarro et al.l(l2010h . 

2.2 Identification of the LMC analog 

We search the Aquarius simulations for accretion events that result 
in systems with kinematics similar to that of the LMC. In particular, 
we look for relatively massive systems (i.e., with masses exceeding 
1% of the main halo mass) that are accreted relatively recently (i.e., 
after z = 1), and that have, at pericenter, distances and velocities 
of order 50 kpc and 400 km/s, respectively. 

Our best candidate is a system that accretes into the main Aq- 
A halo at z = 0.51 (t — 8.6 Gyr, i.e., ~ 5 Gyr before the present 
time). Just before accretion, at Zi^ = 0.9 (tid = 6.6 Gyrs), this 
subhalo has a virial mass of M200 = 3.6 x 10^° Mq. Our analysis 
below focuses on this "LMC-analog" halo and its substructures as- 
sociated to the same friends-of-friends group ("LMCa group", for 
short) in the Aq-A-3 simulation. This level-3 simulation has a par- 
ticle mass nip — 4.9 x lO^M©, and therefo re LMCa is resolved 
with more than 700, 000 particles. SUB FIND jSpringel et al.l200ll) 
is used to identify self-bound substructures within LMCa; at tid 
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Figure 1. Left panel: Time evolution of the galactocentric distance of tlie "LMC-analog" (LMCa) halo (solid black curve) as well as that of its subhalos 
(dashed cui'ves), identified at 4;^ (vertical arrow). Note the complex orbital evolution of associated subhalos, which include a few that are temporarily captured 
by LMCa after being "ejected" from the main halo during the tidal dissociation of a group accreted at t ~ 4 Gyr The LMC analog enters the viiial radius of 
the main halo (marked by the dotted line) at t ~ 9 Gyr and reaches its first pericenter at t ~ 9.6 Gyr After first pericenter the subhalo loses 38% of its mass 
and its orbital apocenter is reduced to ~ 250 kpc, about half of its first turnaround distance. A non-negligible fraction of subhalos (~ 4%) are ejected into 
very eccentric orbits after first pericenter, whereas others settle onto more bound orbits with shorter periods. The main subhalo reaches its second pericenter 
alt ^ 13.3 Gyr still surrounded by some of its most bound companions. Right panel: Radial and tangential velocity of the main subh alo as a function of 
time. Shaded areas correspond to the observed velocity of the LMC (±3o") based on the proper motion measurements of lKallivavalil et aL ! (2006). Blue dotted 
vertical lines indicate the times tip and i2p (near first and second pericenter passages) when the kinematics of LMCa best matches that of the LMC. 



there are more than 250 LMCa subhalos with masses exceeding 
1 X IO^Mq. 

As may be seen from Fig. [T] the LMCa group turns around 
at tta = 5 Gyr {z — 1.3) from a distance of rta = 480 kpc (all 
distances and velocities quoted are physical unless explicitly stated 
otherwise). LMCa reaches its first pericenter at tpcri = 9.5 Gyr, 
with pericentric distance (rpcri ~ 63 kpc) and speed (Vpavi = 355 
km/s) comparable to those of the LMC. The LMCa orbit becomes 
substantially more bound after first pericenter due to dynamical 
friction and tidal mass loss, so that it reaches a distance of only 240 
kpc at its second apocenter, al t — 11.5 Gyr. Its radial period re- 
duced to 3.8 Gyr, it goes through pericenter again at tpcr2 ~ 13.3 
Gyr, with similar pericentric distance and speed as the first (see 
right panel of Fig. [T). 

The best match to the kinematic properties of the LMC occurs 
just after each of these pericentric passages, at times that we will 
denote tip and t2p. These are shown in the right panel of Fig. [T| 
with vertical dotted lines. At these times, the virial mass of the host 
halo is 1.6 and 1.8 x 10^^ Mq, respectively. At t = tip{t2p), the 
satellite is at a distance of 65(69) kpc, with radial velocity V,- — 
78(89) km/s and tangential velocity Vt = 347(302) km/s. In the 
analysis that follows we focus on the properties of LMCa and its 
associated subhalos at tip and t2p. 

These values are in reasonable agreement with the latest LMC 
measurements; in particular, the tangent ial speed at first pericen- 
ter agr ees within la with the results of lKallivavalil et all i200d) 
and of IPiatek et alj j2008h . The velocities are lower by ~ 15% at 
second pericenter, mainly because the apocenter of the orbit has 
been reduced to 230 kpc from the 480 kpc reached at turnaround. 



This argument would appear to favour a first-pericentric interpreta- 
tion for the LMC orbit, in tune with recent suggestions (see, e.g., 
Besla et al. 2007, Boylan-Kolchin et al. 2010, Busha et al. 2010b, 
Tollerud et al. 2011). We urge caution, however, with this inter- 
pretation. The actual pericentric speed will depend sensitively on 
the actual turnaround radius and on the mass of the Galaxy, both 
of which are rather uncertain. Further, we do not expect an exact 
match between LMC and LMCa because of the limited statistics 
that the six Aquarius halos allows. A more robust result seems to 
be that, because of the large pericentric distance, systems analogous 
to the LMC should experience a relatively moderate loss of orbital 
energy, allowing them to return to its second pericenter with similar 
speed to the first. If LMCa had had a slightly larger first-pericenter 
speed then quite possibly the second pericenter would have been 
in better agreement with the LMC than the first. We conclude that 
multiple passages cannot be confidently ruled out by this evidence 
alone. 



3 RESULTS 

3. 1 The Orbit of the LMCa Group 

The left panel of Fig.[l]shows the radial evolution of LMCa and its 
250 most massive associated subhalos (identified at tid within the 
same friends-of-friends group). This shows the complex orbital be- 
haviour of LMCa subhalos as the group gets disrupted in t he tida l 
field of the primary halo. As emphasized by ISales et si. I (l2007ah 
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Figure 2. Aitoff projection of all particles belonging to the LMCa group at two different times, tip (near 1st pericenter, left panel) and t2p (near 2nd pericenter, 
right panel), when the kinematics of the most massive LMCa member (shown by a green circle) matches best that of the Large Magellanic Cloud (see Fig.[T). 
Dots denote all individual particles within the virial radius of the LMCa at the time of identification, tij . Those colored black are still bound to LMCa; grey 
are those that have become unbound as a result of the tidal interaction with the primary halo. Particles in red are those belonging to self-bound substructures 
of LMCa. The positions in the sky of all known Milky Way satellites are also shown a nd labelled in blue. T he orbital path of LMCa is shown by the dashed 
green curve; a magenta curve traces the location of the Magellanic Stream, taken from lNidever et all bOlOl) . The reference frame of the simulation has been 
rotated so that the position in the sky and the proper motion of LMCa match those of the LMC. Note how all particles and substructures associated with LMCa 
stretch along the orbital path. At tip (panel on the left) most associated subhalos lie near the LMC since tidal stripping has not yet had enough time to disrupt 
the system. At t2p (panel on the right) the particles and subhalos associated with LMCa are more widely spread across the sky but still trace the orbital path. 
The LMCa footprint can be used, together with velocity information, to investigate whether other Milky Way satellites have been associated with the LMC in 
the past. See text for further details. 



and lLudlowetalJbOOgI) . many satellites undergo drastic changes 
in energy and angular momentum during pericentric passage, caus- 
ing some systems to achieve escape velocity and to leave the pri- 
mary halo altogether 

Interestingly, the LMCa group itself contains captured "es- 
capees" from another group. As may be seen in the left panel of 
Fig. [T] a number of subhalos belonging to LMCa at tid were ac- 
tually accreted much earlier into the main halo and ejected from it 
when their parent group was disrupted at t ~ 4 Gyr By chance 
these were then temporarily captured by LMCa at tid before plung- 
ing again into the main halo. 

These changes limit the applicability of the usual proce- 
dure of looking for satellite associations by searching for ob- 
jects of common energy and/or angular momentum (see, e.g., 
iLvnden-Bell & Lvnden-Beiilll995t) . Instead of making simple as- 
sumptions we shall use directly the 6D inforination from the simu- 
lations to assess the probability that other Milky Way satellites are 
associated with the LMC. 



3.2 Position and Kinematics of the LMCa Group 

We plot in Fig. [2] the position in the sky (in an Aitoff projection) 
of all particles identified at tid to belong to LMCa. The panel on 
the left corresponds to first pericenter, tip, and that on the right to 
second pericenter, t2p ■ The simulation reference frame has been ro- 
tated so that the position in the sky and the direction of motion of 
the main LMCa subhalo (indicated by a green circle) coincide with 



those of the LMC. The orbital path is shown by a gree n dashed line; 
the li ne in magenta traces the Magellanic Stream dNidever et al.l 
bold) . Particles in black are those currently still bound to LMCa, 
those in grey correspond to tidally stripped LMCa material. Par- 
ticles in red highlight the LMCa substructures that remain self- 
bound. 



Fig.|2]shows the complex footprint on the sky traced by LMCa 
and its debris, and how it changes substantially from first to second 
pericentric passage. At first approach (left panel of Fig. |2), most 
LMCa subhalos cluster around the LMC because tidal stripping 
has not yet had enough time to disrupt the system. Nevertheless, 
there are a few tidally-stripped systems (those more loosely bound 
to LMCa at tid), and they roughly trace the orbit of LMCa along 
a nearly polar circle. The sky is not uniformly populated by the 
stream; for example, because the pericentric passage occurs near 
the Southern Galactic Pole very few LMCa subhalos have a chance 
of reaching the northern Galactic cap by tip. Some systems stray 
from the orbital plane despite their previous association, but these 
cases are rare and correspond to subhalos populating the outskirts 
of the LMCa group at tid- 

The orbital path is still discernible in the tidal stream at second 
pericenter, but the debris populates now both hemispheres in the 
sky. At this time LMCa has lost more than 70% of its original mass 
to tides, compared with only 40% ait = tip. The distribution of 
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Figure 3. Galactocentric distance versus radial velocity for all particles in the LMCa group, compared with those of known Milky Way dwarfs. Color coding is 
the same as in Fig.|2] Left and right panels correspond to t = tip and t = t2p, respectively. For reference, the local escape velocity of a Navarro-Frenk- White 
halo (Navarro et al. 1996, 1997) is indicated by the blue long-dashed curves, for two values of the concentration, c = 10 and 20. Note that the tidal debris from 
LMCa is constrained to a well defined region in the (r, Vr) plane, a fact that may be used to gauge the likelihood of association with the LMC of individual 
dwarfs. 



the debris on the sky would certainly be wider for a progenitor of 
larger mass than the one considered here 

Blue asterisks in Fig. [2] indicate the position of all known 
Milky Way satellites, including the ultra-faint dwarfs discovered 
by SDSS, which are located mainly in the northern Galactic cap 
(see Table |2]for the compilation of values used). Several dwarfs 
lie along the debris path and therefore could potentially have been 
associated with the LMC. These include, att — tip. Sextans (Sex), 
Draco (Dra), Segue 3 (Seg3), Pisces II (PiscII), Sculptor (Scl), For- 
nax (For), Tucana (Tuc), Carina (Car), Leo II, and the SMC. Be- 
cause the debris is more spread out at t2p, the list includes then 
further dwarfs, such as those in the constellation of Leo. 

Of course, true association requires not only coincidence with 
the LMCa tidal stream in projection, but also in distance and veloc- 
ity. We investigate this in Fig.|3] where we plot the Galactocentric 
distance and radial velocity of LMCa particles for t — tip (left) 
and t = t2p (right). Color coding is the same as in Fig.|2] Data for 
all known Milky Way satellites (as listed in Table|2l( are also shown 
in each panel, for comparison. Since the LMCa debris is confined 
to specific regions in the (r, Vr) plane, these data may be used to 
test the association of any individual dwarf with the LMC, subject 
to assuming that the LMC is either on its first or second approach 
to the Galaxy. We explore these associations next. 

3.3 Association with the Clouds 

To be deemed an "LMC associate" a satellite must satisfy at least 
the following three conditions: (i) it should fall in the celestial 
sphere within the footprint of the LMCa group; (ii) it should be at 
a Galactocentric distance consistent with that of LMCa particles at 
the same location in the sky; and (iii) it should have Galactocentric 
radial velocity also consistent with LMCa in the same region. In 
principle a further condition involving the tangential velocity could 



be added to the list, but the preliminary nature of most proper mo- 
tion estimates for Milky Way satellites implies that the results are 
unlikely to be conclusive at this stage (see the Appendix for more 
details). We discuss conditions (i) to (iii) below. 



3.3.1 Sky proximity 

We can quantify the proximity in the sky of any dwarf to the LMCa 
stream by computing the angular distance to the n"' nearest LMCa 
particle in projection and comparing that with the probability of ob- 
taining a similar distance (or smaller) for a random point in the sky. 
Note that the analysis that follows uses all LMCa particles (rather 
than just subhalos) when comparing with the Milky Way satellites, 
since this provides a more complete sampling of the distribution in 
phase space expected for LMC debris. 

The histograms in Fig. |4] show the distribution of the angu- 
lar distance to the 50th nearest LMCa particle, S50, computed for 
10,000 random points in the celestial sphere; the median of the dis- 
tribution is at (5iim = 4° and Sum ~ 2° for the first and second 
pericenter passages, respectively. As may be seen from the inserts 
in Fig. |4] the footprint of the LMCa stream is traced faithfully by 
regions satisfying S50 < <5iim (shown in red) and we shall use this 
condition to decide which dwarfs are likely to be associated with 
the LMC. 

According to Fig.|4] at first pericenter only 8 dwarfs satisfy the 
proximity condition specified above but, because the LMCa foot- 
print covers a much larger fraction of the sky at second pericen- 
ter, 17 dwarfs pass this constraint then. These numbers are reason- 
ably insensitive to our choice of n=50, especially at first pericenter; 
choosing the 100th nearest neighbour results in 9 likely-associated 
dwarfs at tip and 22 at t2p. 
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Figure 4. Distribution function of ^so, the angular size (radius) of tlie circle that contains the 50 nearest LMCa particles in projection, computed for 10, 000 
random points in the celestial sphere. Left and right panels correspond to t = tip and t = t2p, respectively. ^50 is a useful measure of proximity to the LMCa 
stream in the sky for any given point. In the inserted Aitoff maps, red is used to highlight points with S^o < (Sum, where is the median of the random 
distribution. Black points correspond to 550 > <5iim- As demonstrated by the Aitoff inserts, 550 is a useful measure of association with the stream. We show 
in each panel the values of <55o corresponding to each Milky Way dwarf, and retain only those with ^50 < Sn^^ as potentially associated with the LMC. 



3.3.2 Radial velocity 



Figs. |5] and [6] are then used to check which of these candi- 
date dwarfs are also associated with the stream in distance and ve- 
locity. The various panels in these figures show the Galactocentric 
distance and radial velocity of each of the dwarfs in the LMCa sky 
footprint, together with those of all LMCa particles closer than Sum 
in the sky from each dwarf. This enables an intuitive test of which 
satellites have distances and velocities consistent with LMC asso- 
ciation. 

For example, the top left panel of Fig. |5] corresponds to the 
SMC and it shows that there are plenty of LMCa particles with 
distances and radial velocities coincident with the SMC at first 
approach, endorsing a true association between the LMC and the 
SMC. An opposite example is provided by Draco; the data in Fig.|5] 
show that, although it is within the LMCa footprint, all of the par- 
ticles associated with LMCa at Draco's location have discrepant 
velocities and/or distances. This is because Draco is projected onto 
the trailing stream at first pericenter, and one would expect then 
very large negative radial velocities at Draco's Galactocentric dis- 
tance. A similar analysis may be carried out for each dwarf individ- 
ually, but it should be clear from Fig. |5]that, aside from the SMC, 
only Carina and Fornax seem to have a reasonable chance of being 
associated with the Clouds if they are on first approach. 

The situation is less clear-cut at second pericenter, because the 
LMCa footprint is larger, and because the stream now has multiple 
wraps, which allows for a wider range of velocities in a given direc- 
tion in the sky to be consistent with LMCa association. The right 
panel of Fig.|4]indicates that 17 dwarfs pass the sky proximity test; 
these are shown (two per panel except in the middle left) in Fig. [6] 



Inspection of each panel shows that, aside from the SMC, Carina, 
Canes Venatici II, Sculptor, Leo II, Leo IV, and Leo V could in prin- 
ciple be associated with the LMC. The situation of Fornax, Canes 
Venatici I, and Sextans is less clear but one would be hard pressed 
to rule out an association in such cases. Bootes (1 and 2), Segue 3, 
LeoT, Willman I and Coma are all unlikely to be associated with 
the LMC. 



3.3.3 Orbital angular momentum 

The sky proximity and radial velocity constraints discussed above 
are sensitive to the fact that we consider here a single LMC look- 
alike system. This is compounded by the relatively low mass 
of LMCa, compared with the 1 .3 x IO'^^Mq sugges t ed by the 
abundance-matching analysis of iBovlan-Kolchin et alj t201ll) . A 
more massive subhalo may leave a broader footprint on the sky 
and lead to a wider range of radial velocities consistent with LMC- 
associated debris. More conclusive statements about the likelihood 
of association of the candidate satellites identified in the previous 
subsections require accurate measurements of the proper motion in 
order to constrain their tangential velocity and to verify that they 
follow orbits roughly aligned with the orbital plane of the LMC. 
Indeed, as we show below, the direction of the orbital angular mo- 
mentum of a dwarf might be one of the cleanest tests of association 
with the Clouds. 

This is shown in Fig. [T] where we show, in an Aitoff projec- 
tion of Galactocentric coordinates (la, bo), the direction of the or- 
bital angular momentum of LMCa and its associated substructures 
at first and second pericentric passage. Because of the nearly po- 
lar orbit the direction of the angular momentum of most LMCa- 
associated material is roughly on the Galactic plane, pointing in 
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Figure 5. Galactocentric radial velocity and distance for dwarfs (filled squares) deemed possibly associated with the LMC according to the criterion of Fig. |4] 
i.e., (Sso < (Slim at t = tip. Each panel also shows the r and Vr of all LMCa particles (dots) within a circle of radius <5ii„j centered at the position of each 
dwarf. This probes graphically whether the positional association indicated by proximity to the stream in the sky is corroborated by the velocity data. This test 
indicates that the SMC is the only known satellite clearly associated with the LMC if the Clouds are on their first pericentric approach. Aside from the SMC 
only Carina and Fornax seem marginally consistent with an LMC association. On the other hand, this test seems to rale out a possible association for all other 
candidates. 



the direction of the Sun from the Galactic center; i.e., la = 180°, 
be ~ 0°. This tight alignment is preserved at second approach, 
albeit with larger scatter. 

In Galactocentric Cartesian coordinates, with the X-axis 
pointing away from the Sun, Y-axis defined positive in the direc- 
tion of Galactic rotation, and positive Z-axis in the direction of 
the Galactic North Pole, this implies that the X-component of the 
orbital angular momentum (jx) of associated satellites should be 
negative and much larger in magnitude than jy or jz- This may 
be seen in Table [T] which lists the components of the unit vec- 
tor identifying the direction of the (average) angular momentum of 
particles associated in the sky with the candidate dwarfs identified 
in the previous subsection. Note that, with no exception, the angu- 
lar momentum points clearly toward —X, the anti-Galactic center 
direction. 

This result makes strong predictions regarding the tangential 
velocity of the candidate satellites, which can be checked against 
observation for the few satellites with available proper motions 
(SMC: Kallivayalil (2006b), Carina: IPiatek et alj j2003l) , Fornax: 



iPiatek et al. l and Sculptor: IPiatek et al.l ( |2006|) ). Inspection 

of Table [T] and Fig. |7] shows that, of the four satellites with pub- 
lished spatial velocities, only the SMC appears associated with the 
LMC. None of the other three (Carina, Fornax, and Sculptor) seems 
obviously associated with the Clouds according to this test (see last 
three columns of Table |2]l. In hindsight this not entirely surpris- 
ing, given the very dissimilar chemical enrichment patterns and gas 
content of the LMC and SMC compared with other Galactic satel- 
lites (Mateo 1998, CaiTera et al. 2008, Harris et al. 2009, Kirby et 
al. 201 la,b). Explaining what drives the diversity in star formation 
history, metal enrichment, and gas fractions of Galactic satellites 
remains a prime challenge for dwarf galaxy formation models . We 
hasten to add, however, that, as the recent revision to the proper mo- 
tion of the LMC illustrates ( Kallivayalil et al. 2006), proper motion 
measurements are exceedingly difficult, and hence our conclusion 
should be revisited when new, more accurate data become avail- 
able. 
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Figure 6. Same as Fig.|5]but for the second pericenter passage. Because the stream covers a wider area in the sky, the number of candidate dwarfs increases. 
We show two dwarfs per panel; arrows indicate systems that lie beyond the plotted region. In each panel red (blue) dots are stream particles in the circle 
around the dwarf indicated by a red (blue) square. This exercise confirms the association of the SMC, and suggests that several other dwarfs have kinematics 
consistent with association with the LMC. See text for further discussion. 



4 SUMMARY AND CONCLUSIONS 

We use cosmological N-body simulations from the Aquarius 
Project to study the orbit of the LMC and its possible association 
with the SM C and other Milky Way satellites in light of new proper 
motion data jKallivavaUl et alj200d . |Piatek et alj2008t) . We search 
the simulations for LMC dynamical analogs; i.e., accreted subha- 
los with pericentric distance (~ 50 kpc) and velocity (~ 400 km/s) 
matching those of the Clouds. 

One suitable candidate (LMCa) is a 3.6 x 10^°Mq system 
accreted at2;'^0.5(t = 8.7 Gyr) by Aq-A, a halo that, at 2; = 
has a virial mass of 1.8 x lO^^M©. LMCa turns around from a 
distance of 480 kpc at tta ~ 5 Gyr (z = 1.3), accretes into the 
Milky Way halo at t = 8.6 Gyr (z = 0.5), and completes two 
pericentric passages by 2; = 0. 

We use the positions and velocities of particles belonging to 
LMCa before infall in order to trace the orbital evolution of LMC- 
associated satellites and to inform the analysis of the likelihood 
that other Milky Way satellites were accreted in association with 
the Magellanic Clouds. Our main conclusions may be summarized 
as follows. 

Near each pericentric passage the kinematic properties of 



LMCa match approximately those of the LMC. This implies that 
(i) the orbit of the LMC is not particularly unusual given the halo 
virial mass, and that (ii) it is difficult to decide, using only kine- 
matical data, whether the LMC is on first approach or has already 
completed a full orbit. 

If the LMC is on first approach, then most of its associated 
subhalos should be tightly clustered around its location. Although 
rare, some LMC-associated systems may still be found well away 
from the LMC but along the orbital path of the group. Since none of 
them has completed a single orbit there are strong position- radial 
velocity correlations that may be used to identify which satellites 
might have been accreted together with the LMC. 

Of the known Milky Way satellites only the SMC is clearly 
associated with the LMC. A case can also be made for Fornax, Ca- 
rina, and Sculptor, but it is not a particularly compelling one. This 
is specially true when considering available proper motion data, 
which suggest that the orbital planes of these three satellites are not 
aligned with that of the Clouds. 

If the LMC is near its second pericenter then several further 
dwarfs qualify for association. Leo II, Leo IV and Leo V, in par- 
ticular, show strong spatial and velocity coincidence with the tidal 
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debris from LMCa, making them prime candidates for past associ- 
ation witii tlie LMC. Persuasive, but hardly conclusive, cases can 
also be made for a handful of other dwarfs, such as Canes Ve- 
natici II and Leo I. These tentative associations may be firmed up 
or refuted using the full spatial velocity, for which our simulations 
make a strong prediction: it must be such that the direction of the 
orbital angular momentum should point unambiguously in the anti- 
Galactic center direction. 

We expect then few, if any, known Milky Way satellites to 
be associated with the Clouds, especially if they are on first ap- 
proach. The simulations, however, make a very clear prediction in 
that case: most satellites associated with the Clouds have yet to dis- 
perse and should therefore be very near them. How many could 
we expect? The number depends strongly on the virial mass of the 
Clouds and is, therefore, highly uncertain. If the mass is as high 
as 1.3 X IO^Mq, as suggested by the abundance-matching anal- 
ysis of iBovlan-Kolchin et alj <201lh . then we would expect it to 
have of order 7 subhalos w ith peak circular velocities exceeding 20 



km/s I Springel et al.ll2008h . This, according to the recent model of 
iFont et alj ( l201lh . is the minimum halo potential depth required to 



host a luminous dwarf. 

Although the numbers seem modest, one may be encouraged 
by the fact that the LMC does have at least one companion (the 
SMC), so the possibility that they are part of a larger "Magel- 
lanic Galaxy" should not be dismissed too quickly. Surveys de- 
signed to target the sky around the LMC (SkyMappe fl, MAPS 



signed to target the sky around the LMC (akyMappe p MAF^ 
jNidever et al.l l20Ilh . NOAO Outer Limits Survey dSaha et al.l 
l20IOh ) should help to unravel the history of the LMC and its com- 
panions. The surroundings of the Clouds might be hiding a trove of 
new Milky Way satellites awaiting discovery. 
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Figure 7. Direction of the orbital angular momentum of LMCa at first and second pericenter passages, compared with those of dwaifs with pubHshed proper 
motions (SMC, Fornax, Sculptor and Carina). The directions are shown in an Aitoff projection of Galactocentric coordinates {Iq and be)- The color coding 
is the same as in Figs.[2]and [5] with the green cross indicating the orbital pole of the central subhalo in LMCa. Note that the viewing perspective is different 
from that in Fig.|2] I > 180° in the fornier coiTcspond to negative longitudes here. All satellites associated with the LMC before infall are expected to have 
their orbital poles in a well defined region of the sky centered approximately at {Iq, 6g) ~ (170° , — 10° ). Of the 4 dwarfs with measured full spatial motions 
only the SMC seems obviously associated with the LMC. 



Table 1. Cartesian components of the unit vector that characterizes the direction of the (average) angular momentum of LMCa particles near each candidate 
LMC-associated satellite, according to the discussion of Fig.|3]and Fig.|6] For each dwarf, first and second rows Hst values ait = tip and t = t2p , respectively. 
Note that, because of the nearly polar orbit of LMCa, the angular momentum points in all cases in the —X direction, i.e., to the Sun from the Galactic center. 
The angular momentum direction is also listed for satellites with published proper motions (third row), and may be used to assess their possible association 
with the Clouds. 
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Table 2. List of propeities of the Milky Way satellites used in Figs.|2]and [3] The last three columns provide information about the likelihood of association 
of each dwarf with the LMC according to the sky proximity (Sec.[3]TT), radial velocity (Sec. |3.3!2) and angular momentum (Sec. |3.33l criteria. Pairs of 
symbols show the results for first and second pericenter (in that order). A "check-mark" (/) indicates good or moderate agreement between the properties of 
a given dwarf and expectations from the simulated LMC analog. If, instead, there is a clear mismatch, the table displays "x". Question marks (?) in the last 
column denote the lack of available proper motions to compute the orbital angular momentum of these dwarfs. Previous association with the LMC requires 
the simultaneous fulfillment of the three conditions (see text for more detail). Note that the tests ai'e applied recursively: the sky proximity is computed for all 
dwarf galaxies, but only those with a positive match are considered for the radial velocity and angular momentum tests. Therefore, after a " X " symbol a d warf 
is removed from the list of possible LMC companions, and a minus si gn "— " fil ls the remaining columns. References are as follow s: ri1:lMate o'l'l998'), [2]: 
[simon & Geha (20()7a), [3]: Belokurov et al. (2010), [4]: Belokurov e t ai] j2007l) [5]:|^lokurov et al. (2009), [6]: Simon & Geha ( l2007tl) . [71: Mart in et all 
120071) , ISL IKochetaU <2009l) , r91: IIrwin et alH2007l) . llOk lWillman et alj )2005h , fl lli lfielokurov et ali <2008i) . 
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APPENDIX A: MAGELLANIC GALAXK-DEBRIS IN THE 
LMC FRAME 

Eor completeness. Fig. lAll shows the distance-velocity plane of 
LMCa debris in a coordinate system centered at the LMC. As the 
system orbits within the host potential, mass is lost to tides and the 
material initially associated to the LMC group gets progresively 
unbound. Substructures that remain bound to the LMC are shown 
in black, unbound material in grey. LMCa substructures are shown 
in red. The SMC falls in the "unbound" region at either first or 
second pericenter, which suggests that the LMC-SMC might not 
longer be a bound pair. This conclusion, however, is sensitive to 
the mass of LMCa which, as discussed in the text, is likely smaller 
than the true LMC mass. Once reliable three-dimensional velocities 
become available for more dwarf galaxies, their distribution in the 
LMC-centered phase space may be used to place further constraints 
on their association to the Clouds. 
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Figure Al. Total 3D velocity, V^j^"^ , versus distance r^'^''^' , where upperscripts indicate that they are computed with respect to the LMC center Color 
coding is the same as in the previous figure. The green asterisk shows the position of the SMC in this plane according to data from Kallivayalil et al. 2006b. 
The transition from black (bound) to gray (unbound) dots can be used to infer the instantaneous escape velocity of the LMCa system. For comparison, blue 
curves show Vesc for an NFW halo with the mass of LMCa at the time of infall, M200 = 3.6 X IO^'^Mq, assuming two different concentrations c=10,20. 
The effect of tides due to the host potential can be seen from the comparison between these curves and the velocity of the bound (black) particles; in particular, 
the mass loss experienced between first and second approach is reflected by the smaller area (and lower velocities) covered by black particles in the right panel 
(second peiicenter) compared to the left (first pericenter) . 
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